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Majorana quasiparticles localized in vortex cores of a chiral p-wave superconductor hybridize with 
one another to form bands in a vortex lattice. We begin by solving a fully microscopic theory 
describing all quasiparticle bands in a chiral p-wave superconductor in magnetic field, then use 
this solution to build localized Wannier wavefunctions corresponding to Majorana quasiparticles. A 
low-energy tight-binding theory describing the intervortex hopping of these is then derived, and its 
topological properties—which depend crucially on the signs of the imaginary intervortex hopping 
parameters—are studied. We show that the energy gap between the Majorana bands may be either 
topologically trivial or nontrivial, depending on whether the Chern number contributions from the 
Majorana bands and those from the background superconducting condensate add constructively or 
destructively. This topology directly affects the temperature-dependent thermal Hall conductivity, 
which we also calculate. 

PACS numbers: 


It has long been known that an isolated vortex core 
in a chiral p-wave superconductor can host zero-energy 
Majorana modes m- Such topological superconductiv¬ 
ity may be present in Sr 2 Ru 04 El IS], and also has ana¬ 
logues in the Moore-Read fractional quatum Hall state 
mn] and in the A-phase of superfluid ^He [5|. When 
multiple vortices are brought close together, the Majo¬ 
rana zero modes hybridize, leading to states at positive 
and negative energy. From studying pairs of vortices, it 
has been shown that the sign as well of the amplitude of 
this energy splitting depends on the intervortex distance 

[3iini. 

Just as there are two possible signs of the energy split¬ 
ting for a pair of vortices, there are two topologically 
distinct possibilities for the energy gap between the Ma¬ 
jorana bands in a system consisting of many vortices m 
[Bj. In this work we provide a way in which to compute 
the topological properties of the ground state starting 
from a microscopic Hamiltonian describing paired elec¬ 
trons on a square lattice. We show how the topologi¬ 
cally trivial and nontrivial states can be understood as 
arising due to destructive or constructive addition of the 
Chern number contributions from the Majorana bands 
and from the background superconducting condensate, 
respectively. We further show that these different topo¬ 
logical states can lead to different behaviors of the intrin¬ 
sic contribution to the thermal Hall conductivity, which 
we compute explicitly from the microscopic model. 

Several authors have recently studied periodic arrays 
of vortices starting from a continuum model of a chi¬ 
ral p-wave superconductor [BHB] . Unlike Ref. [T3], we 
do not find any zero-energy flat band with zero Chern 
number. This disagreement arises due to the neglect of 
the magnetic field in that work, as well as the different 
choice of the spatial profile of the complex phase of the 


superconducting order parameter. Our work differs from 
Ref. [16] in that we do not find any strong anisotropy 
emerging in our low-energy theory describing Majorana 
quasiparticles. Further, the authors of that work find 
a single edge mode in their low-energy theory, whereas 
an even number always emerges from our theory. At 
present, we do not understand the precise reason for this 
disagreement. Our results extend those of the semiclassi- 
cal calculation in Ref. m by providing a fully quantum 
treatment, which is necessary in order to address effects 
associated with Berry phases and topological properties. 

Finally, very recently another work that has some over¬ 
lap with our results has appeared m- As in our work, 
the authors begin from a lattice Hamiltonian describing a 
chiral p-wave superconductor in magnetic field and solve 
for the full quasiparticle bandstructure. They then ob¬ 
tain a tight-binding description of the Majorana bands by 
fitting the hopping parameters to this full solution. This 
is a different procedure from our work, in which deriving 
a tight-binding model from localized Wannier functions 
allows us to determine the signs of these hopping param¬ 
eters rather than just the absolute values. As we will 
show below, these signs are what determine the topolog¬ 
ical properties of the Majorana bands, and hence have 
a direct affect on the Berry curvature and thermal Hall 
conductivity. The calculation of these quantities is the 
main result of our work and did not appear in Ref. m 


MICROSCOPIC LATTICE MODEL 

Following the approach of Refs. mills], we use the 
following Hamiltonian to describe a superconductor on a 
square lattice in a magnetic field: 
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^r,r+5Cp^^Cr+5,(T ~t” ^r,r+(5 


Here Cr^a- is the electron operator, with implicit sum¬ 
mation over spin a i- The hopping, which in¬ 
cludes a Peierls phase factor due to the magnetic field, is 
given by t^^r+s = where = Try^/(l>o, 

^r.r-i-y = — 4>Q = hcje is the elementary flux 
quantum, and $ is the magnetic flux through each pla- 
quette. The superconducting pairing term is Ar^r-i-a = 

exp ■ Tor chiral p-wave pair¬ 

ing, the amplitudes are given by (Ar,r±x, Ar.r±y) = 
(±A,±iA) or (±A,=FiA) for ^^Px+iPy' or ^^Px—iPy^ pair¬ 
ing, respectively. As discussed in previous works usms], 
the phase of the superconducting order parameter 9{r) is 
chosen to be a solution to the continuum London equa¬ 
tions: 


V X V0(r) = 27r.2^ (5(r — r^-), 
j 

V26»(r) = 0. 


( 2 ) 


Finally, the chemical potential term in general includes 
a spin-dependent Zeeman term: = p zL hz- When 

written in the Bogoliubov-de Gennes language, the Zee- 
man term is proportional to the identity matrix in the 
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Hamiltonian (see ([^ below), so that hz merely shifts 
the energy bands as a chemical potential would for an or¬ 
dinary fermionic system without superconductivity. We 
consider the effects of this shift when calculating the ther¬ 
mal Hall conductivity below. 

We consider a magnetic unit cell of Lx x Ly sites, with 
each unit cell containing two vortices, as shown in Figure 
The superconductor is assumed to be strongly type- 
H, with a sufficiently small vortex core size and suffiently 
large penetration depth that the amplitudes of both the 
order parameter and magnetic field are constant in space. 
By employing a singular gauge transformation USHS], 
the Hamiltonian Q can be made periodic with spatial 
period Lx(y) along the x (y) direction, thereby enabling 
the use of Bloch’s theorem. Performing the following 
particle-hole transformation: 







( 3 ) 


where kx,y € new fermionic oper- 

X XI ^y 

ators satisfy ^r+R,o-(k) = ^r,CT(k), the Hamiltonian (IB 
becomes 


'^ = E E( E -f iL.C. - (k) V’r,<T (k) 

r^u.c. k \5—x,y 

+ E E E [Aiy5(e*‘^'V^,t(k)V’r+u(k) + e-*'^'VET(k)V’rq(k))+H.C.' 


r^u.c. k S—x,y 


( 4 ) 


k r,r'^u.c. 


where = tK,r+s = = 

and = zLp -|- hz- The 

branch cuts in are chosen to connect vortices pair¬ 

wise within a unit cell, as shown in Figure (Further 
details regarding the treatment of branch cuts can be 
found in Ref. [H].) 

Figure shows the band structure obtained by diag¬ 
onalizing the Hamiltonian The spectrum features 
a large gap with energy ~ A, with two bands near 
e = 0 arising from the intervortex tunneling of Majo- 


rana zero modes. The bandwidth of these low-energy 
bands decreases exponentially with the distance between 
vortices. Upon closer inspection, one finds that the Ma- 
jorana bands feature two gapped Dirac cones, which is 
expected from the fact that the imaginary hopping pa¬ 
rameters describe a tt-Aux model, as pointed out previ¬ 
ously [22] and described in detail below. 

Using the numerically determined eigenstates of the 
Hamiltonian Q, the integrated Berry curvature summed 
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FIG. 1: Square vortex lattice, with arrows indicating the sign 
of the imaginary nearest-neighbor hopping of fermionic quasi¬ 
particles between vortices. The dotted lines indicate branch 
cuts connecting the two vortices within each magnetic unit 
cell. 


over all occupied bands is given by [HIISIIM] 
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In this equation, m and n denote the quasiparticle bands, 
and the summation is over all m and n subject to the 
constraint Em,{k) < ^ < En{k). Methods for evaluat¬ 
ing this expression efficiently were discussed in Ref. m- 
Figure]^ shows the behavior of ^xy{i)- Near the top and 
bottom portions of the energy spectrum, the quasiparti¬ 
cle bands are either electron- or hole-like, so that cF^yiC) 
just sums up the Chern numbers ±1 as ^ passes through 
each quasiparticle band. A large jump from positive to 
negative values occurs at the van Hove singularity, where 
there is a pile-up of Berry curvature. For |^| < A, dxyiC) 
is suppressed due to the fact that the quasiparticle states 
are superpositions of electron and hole states, and hence 
carry little or no Berry curvature. Closer inspection, 
though, reveals an interesting structure in this region. 
Rather than completely vanishing, the integrated Berry 
curvature within the large energy gap takes a value ±I 
for px ± ipy pairing, reflecting the topological nature of 
the superconducting state. At the lowest energies, how¬ 
ever, the Chern number becomes either 0 or ±2. As we 
shall show below, this reflects the fact that the Majo- 
rana bands carry Chern number ±1. One thus sees that 
the topological nature of the gap between the Majorana 
bands depends on the interplay of the topologies of the 
Majorana bands and of the background condensate. 


FIG. 2: (a) Quasiparticle energy bands of a p-wave super¬ 
conductor with Px — ipy pairing in a perpendicular magnetic 
field, for a square vortex lattice with intervortex distance 
Lx = Ly = 12, pairing strength A/t = 0.5, chemical potential 
p/t = —2, and Zeeman shift hz = 0. (b) Detailed view of the 
two Majorana bands, where the dashed lines are obtained by 
solving the full Hamiltonian Q, while the (virtually indis¬ 
tinguishable) solid lines come from the effective tight-binding 
theory with first- and second-neighbor hoppings, (c) Mag¬ 
netic Brilluoin zone, showing the path along which the en¬ 
ergy bands are plotted, (d) The Majorana bands feature two 
Dirac cones along the line ky = —kx, which are gapped due 
to next-nearest neighbor hopping. 


EFFECTIVE TIGHT-BINDING MODEL FOR 
MAJORANA FERMIONS 


In order to better understand the topological nature of 
the Majorana bands, we proceed to derive a 2 x 2 tight- 
binding Hamiltonian describing the hopping of Majorana 
quasiparticles between vortices. Given Bloch eigenstates 
'0nk(r, cr) of the Hamiltonian (|^, Wannier functions can 
be defined as 


Wa(r-R-rQ,,cr) 
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The coefficients appearing in are thus given by 


= (V'nkl^ak) = ^ Ck (j*, ct)V^ ak (r, Cr) . (10) 

rtr 


With the localized Wannier functions in hand, we are 
now in a position to derive a low-energy tight-binding 
Hamiltonian to describe the hopping of quasiparticles be¬ 
tween vortex sites. We begin by expanding the fermion 
operators appearing in Q (putting a hat on the opera¬ 
tors to avoid confusion) in a basis of Bloch eigenstates: 

tpra = e“*'‘''’Vir.<7(q) = V’kn(r, Cr)a„(k). (11) 

kn 

We can then invert to obtain the Bloch states in terms 
of Wannier functions: 


FIG. 3: The integrated Berry curvature, summed over bands 
up to energy with same parameters as in Figure The 
lower panel provides a detailed view near ^ = 0, showing that 
the Chern number within the large energy gap is u = — 1, 
while the Majorana bands carry Chern numbers u = ±1. 


V^„k(r, a) = - R - r„, a),(12) 

R,q 

and rewrite the Hamiltonian as 


where n is a quasiparticle band index, TVuj. is the number 
of points in momentum space, and a = A, B denotes the 
two vortex sublattices. The vectors R and Tq, point to 
a magnetic unit cell and to the position of a vortex of 
type a within the unit cell, respectively. In building a 
low-energy theory, we keep only the values of n corre- 

(k) 

sponding to the two Majorana bands. In this case Uan 
is an arbitrary 2x2 unitary matrix at each k, which we 
want to choose such that each Wannier function is local¬ 
ized at one of the vortex cores located at r = R-f r^. This 
can be accomplished straightforwardly using projection 
operators starting from a trial wavefunction |251l26j . The 
first step is to define initial guesses I^rq) for the Wannier 
functions, which we take to be real Gaussian functions 
localized near each vortex core. We then project these 
trial wavefunctions onto the Bloch functions for the two 
Majorana bands: 

l^ak) — ^ ^ |'0mk) (V^mk|^/R—0,a)- (7) 


EE 

r<7 r'a' 

= EEEE“"W“"'(^') 

r<7 v' g' nk n'k 

X V'nk(i', o)U(va, r'cr')V'n'k'(i’', cr'). 


where a') = 'Hk=o(rcr, r'cr') is the Hamiltonian 

in real space before applying Bloch’s theorem. Using 


(121 and defining the new fermion operators that create 


quasiparticles on vortex sublattice a 


dl(k) = 5]t/WaUk), 


(14) 


the Hamiltonian becomes 


^-EE fl)j(k)'H„a'(k)da'(k), (15) 

OiOc' k 


where 


Unlike the original Bloch functions, which have an ar¬ 
bitrary phase at each k, these states will have phases 
that are smooth as a function of k. The next step is to 
orthonormalize these states to obtain pseudo-Bloch func¬ 
tions: 

I^Qk) = ^ |<?i'Qk)(5'k ^^^)na, (8) 

n 

where {Sk)na = Ju,c.i^rik\4>ak) (the integral refers to av¬ 
eraging over the unit cell). The Wannier functions are 
then given by 

k 


y-aa' (k) — E E E ' 


,ik- (AR+r^/ — Fq, ) 


ra y'g' AR ( 16 ) 

X w* (r — r„, u)'R(va, v' G')wa'ij' — AR — Tq,/, ct'). 


Due to the localization of the Wannier functions, only a 
limited number of terms with small AR will give signifi¬ 


cant contributions to (16). 


While the sums in (16) could be computed using 
the Wannier functions constructed in (|^, a simpler ap¬ 
proach is to reexpress these Wannier functions in terms 
of the Bloch function eigenstates of Hk(rcr, rV') = 
g-ik r-^(rij, r'(j')e*k r ^ obtains 
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rcr AR qn q'n' 

— ^ gik (AR+r„,-r„) ^ ^-iq-(AR+r„/-r 
AR qn 


g-- ..,^^-,q.(AR+r„,-r„)(y(q))*^(q)^^^^ 


(17) 


Given the Um determined by the Wannier functions in 


(10 1 and the energies for the Majorana bands shown 


in Figure this sum can be easily evaluated numeri¬ 
cally to determine the tight-binding parameters describ¬ 
ing hopping between vortex cores separated by distance 
AR -I- Ta' - Ta- 

For a square vortex lattice, the nearest- and next- 
nearest-neighbor hoppings are all imaginary, with signs 
of the nearest-neighbor hoppings as shown in Figure 
This corresponds to a tt-Aux state, in which each pla- 
quette is penetrated by a half quantum of magnetic flux 
[niiTHin]- The low-energy effective Hamiltonian takes 
the form 


TiMik) =-to sin 


CTl 


+ to cos cr 2 

-I-ti[sin(fc 2 ,La;) - sin(fcyLy)]CT 3 . 


(18) 


The real parameters to and ti correspond to nearest- and 
next-nearest-neighbor hopping, respectively [32] • Diago¬ 


nalizing (181 leads to energy bands with dispersion 


4 = ± 


tn COS^ 


-htgsin^ 


tl[sin{ka:Lx:) - sm{kyLy)Y 


1/2 


(19) 


which describes two massless Dirac cones at the points 
{kx,ky) = ±( 7 r/La;, — 7 r/Ly) when ti = 0 , with these 
cones becoming gapped when ti 7 ^ 0, as shown in Figure 

m [ssj- 

As expected from previous studies laiiQKiT], the tight- 
binding parameters to and ti oscillate as a function of 
kpL, where kp the Fermi wavevector and L = 
is the magnetic length. Because kp is determined by the 
chemical potential, we show the oscillations of to and ti 
as a function of /r in Figure We also find, as shown 
in Figure l^b), that the ground state energy (defined as 
the sum of all quasiparticle energies over occupied bands) 
becomes lower for px — ipy pairing than for px + ipy as 
/i is increased when the magnetic field is along the +z 
direction (as determined by the Peierls phase factors in¬ 
troduced below 0 ). This illustrates that, although the 
two possible chiralities of the order parameter couple to 
the magnetic field with opposite sign, determining which 
pairing state is energetically favorable is rather subtle. 




p/t p/t 



pit 

FIG. 4: (a) Oscillations in the nearest- and next-nearest- 
neighbor hopping parameters for the Majorana bands de¬ 
scribed by the low-energy effective Hamiltonian ( |18[ ). Solid 
(dashed) lines correspond to to (ti, multiplied x lO), and blue 
(red) lines correspond to Px — iPy {Px + iPy) pairing, with 
magnetic length Lx = Ly = 16, chemical potential p/t = —2, 
and pairing strength A/t = 0.5. (b) Energy difference per site 
between the two pairing states, showing that px — ipy pairing 
is favored for larger values of p. (c) Phase diagram showing 
the total Chern number in the gap between Majorana bands, 
which is determined by the type of pairing and the sign of ti. 
The pairing is Px + ipy to the left of the red line and Px — ipy 
to the right. 


with the result depending on nonuniversal parameters 
that influence the detailed structure of the quasiparticle 
bands. (Note, however, that we have not determined the 
value of A from a fully self-consistent calculation, which 
should be performed in order to make a more definitive 
comparison between different candidate order parame¬ 
ters.) 

Writing the Hamiltonian as ’H(k) = ^i(k)cri, one 

finds that the Majorana bands have nontrivial Berry cur- 
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vature [50] : 


rik 


1 



5h(k) ah(k) 

dk^ dky 


( 20 ) 


with the upper (lower) sign corresponding to the higher- 
(lower-) energy band. Integrating the Berry curvature 
over momentum gives the Chern number of each Majo- 
rana band: 


V = =Fsgn(ti), 


( 21 ) 



where again the upper (lower) sign corresponds to the 
higher- (lower-) energy band. 

Finally, although the focus so far has been on the 
square vortex lattice, we note that the effective Hamil¬ 
tonian (18) can also describe non-square vortex lattices 
with Lx ^ Ly if the next-nearest-neighbor hopping terms 
are allowed to be different along different directions, 
i.e. the 173 term is replaced by ti sixi(kxLx) — t'-^ s\Ti{kyLy). 
Such anisotropic vortex lattices do not in general have the 
Dirac cones that appear in the square lattice, even when 
further-neighbor hopping is tuned to zero. The Majorana 
bands still have Chern number ±1. In the special case 
of a triangular vortex lattice [34] . where LxjLy = VS, 
one has (ti,ti) —t (0,to); and the Chern numbers of the 
Majorana bands are simply determined by the sign of tg, 
as shown in previous work m- 


EDGE STATES AND THERMAL 
CONDUCTIVITY 

When the topologically nontrivial Majorana bands are 
considered within the background of the chiral p-wave 
condensate, two possibilities arise: either the Chern num¬ 
bers of the Majorana bands add constructively with those 
of the condensate, giving overall Chern number = ±2, 
or they add destructively, giving v = 0. A phase diagram 
showing these cases is shown in Figure |^c). The edge 
states for the two cases are shown in Figure]^ from which 
it can be seen that the number of edge modes corresponds 
to the Chern number in each case. In particular, in both 
cases there is a single chiral edge mode in the large energy 
gaps above and below the Majorana bands, reflecting the 
chiral nature of the superconducting condensate. (These 
gaps could be accessed by shifting the Majorana bands 
with a Zeeman term, which is not included in the figure.) 
Figurej^a) shows a case in which the Chern number aris¬ 
ing from the Majorana bands cancels with that of the 
background condensate, so that there are no edge modes 
near energy £ = 0. In contrast. Figure [^b) shows a case 
in which there are two chiral edge modes near £ = 0 , with 
one arising from the superconducting condensate and the 
other arising from the topological Majorana bands. 

Although these edge modes do not contribute to elec¬ 
trical transport due to the fact that they occur within 


0.2 

0.1 

0.0 - 
t 

-0.1 


0 

iO 

FIG. 5: Energy eigenvalues for a chiral p-wave superconduc¬ 
tor (with px — ipy pairing and magnetic length = Ly = 8) 
in the vortex state on a cylinder with periodic boundary con¬ 
ditions along the j/-direction and open boundary conditions 
along the a:-direction. The right panels show the correspond¬ 
ing integrated Berry curvature summed over occupied bands, 
(a) With p/t = —2, the gap between Majorana bands is topo¬ 
logically trivial due to cancellation of the Berry curvature of 
the Majorana bands with that of the superconducting conden¬ 
sate, so that no edge modes are present. At energies above 
and below the Majorana bands, however, states localized at 
the right (red) and left (blue) edges are present due to the 
nontrivial topology of the background superconducting con¬ 
densate. (b) With p/t = —1.5, the gap between Majorana 
bands is topologically nontrivial due to constructive addition 
of the Berry curvature of the Majorana bands with that of 
the condensate, so that two chiral edge modes are present . 



a superconducting state, they do contribute a thermal 
current in the presence of a temperature gradient. The 
thermal Hall conductivity is given by [Biiaisi] 




dnF^\ 

j 


CTxy(0, (22) 


where npiO = l/(e^/^ -I- 1) is the Fermi occupation fac¬ 
tor. The thermal Hall conductivity is thus determined by 
the integrated Berry curvature, convolved with a ther¬ 
mal factor that is peaked near ^ = 0. This quantity 
is plotted in Figure In Figure [^a) we show the re¬ 
sult without Zeeman shift (hz = 0), and in this case 
Kxy{T)/T is determined by the d'xy(0 near energy ^ = 0 
at low temperature. At higher temperature, the ther¬ 


mal factor in (22) broadens beyond the width of the Ma¬ 


jorana bands, so that Kxy(T)/T is instead determined 
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FIG. 6: Thermal Hall conductivity divided by temperature 
for A/t = 0.5, Lx,y = 8, and Px — ipy pairing. The black 
dotted line corresponds to the quantized value for a system 
with Chern number v = —1. Solid red and dashed blue lines 
show the results for the cases where the gap around energy 
^ = 0 has Chern number v = —2 (at p/t = —2) and i/ — 0 
(at p/t = —1.5), respectively (see phase diagram in Figure 
gc))- Insets show the corresponding axy{^), which deter¬ 
mine Kxy(T) via ( |22[ ). (a) With no Zeeman coupling, the low- 
temperature behavior is determined by a^yiC) near ^ ~ 0. 
(b) With Zeeman coupling hz = O.lf, the low-temperature 
behavior is determined instead by OxyiX) near ^ ^ hz- 


by ^xyiO energies above and below the Majorana 
bands. Figure |^b) shows Kxy{T)/T with Zeeman shift. 
As noted above, the effect of Zeeman coupling is sim¬ 
ply to shift the overall energy of the quasiparticle bands, 
so that d'xyi^) —>■ ^xy{i — hz), and KxyiT) is deter¬ 
mined via (22) by the value of dxy near ^ ^ hz- In 
this example, we choose the Zeeman energy hz to be 
half of the cyclotron energy, as it is for free electrons: 
hz = hoJc/‘2 = 2'Kt/[LxLy), where we have expanded the 
tight-binding dispersion near its minimum to obtain the 
effective mass m = h?/{2t)- In this case Kxy{T)/T is 
determined by axyii ~ hz) at low temperatures, while 
it shows a modest change as the thermal factor in (22) 
broadens at higher temperatures, increasing or decreas¬ 
ing to reflect the behavior of hxyi^) near ^ ~ 0. 


CONCLUSION 

The results of our work may have relevance for the can¬ 
didate topological superconductor Sr 2 Ru 04 , for which 
very clean samples are available and thermal conductiv¬ 
ity measurements are feasible. Because it is a multiband 
superconductor, however, a more realistic lattice model 
should be used if quantitative comparison is to be made. 
Our results may also be relevant for the Moore-Read frac¬ 
tional quantum Hall state, which can be thought of as a 
gas of composite fermions with chiral p-wave pairing. 

In conclusion, we have shown that the topological na¬ 
ture of the vortex state of a chiral p-wave superconduc¬ 
tor at low energies is far from being obvious, generally 
requiring a complete microscopic calculation to deter¬ 
mine reliably. In particular, the result depends on the 
constructive or destructive addition of Chern numbers 
associated with the Majorana bands and the supercon¬ 
ducting condensate. More broadly, our work illustrates a 
way in which thermal conductivity can be a useful probe 
for obtaining information about the quantum wavefunc- 
tion describing an interacting many-body system. Such 
a probe is especially useful in a superconductor, where 
electrical transport measurements may not yield useful 
information. 
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